October 2012 


LOUGHNER ET AL. 


1775 


Roles of Urban Tree Canopy and Buildings in Urban Heat Island 
Effects: Parameterization and Preliminary Results 

Christopher P. Loughner 

Earth System Science Interdisciplinary Center, University of Maryland, College Park, College Park, 
and NASA Goddard Space Flight Center, Greenbelt, Maryland 

Dale J. Allen and Da-Lin Zhang 

Department of Atmospheric and Oceanic Science, University of Maryland, College Park, College Park, Maryland 

Kenneth E. Pickering 

NASA Goddard Space Flight Center, Greenbelt, Maryland 

Russell R. Dickerson 

Department of Atmospheric and Oceanic Science, University of Maryland, College Park, College Park, Maryland 

Laura Landry 

Maryland Department of the Environment, Baltimore, Maryland 


(Manuscript received 28 October 2011, in final form 16 April 2012) 

ABSTRACT 

Urban heat island (UHI) effects can strengthen heat waves and air pollution episodes. In this study, the 
dampening impact of urban trees on the UHI during an extreme heat wave in the Washington, D.C., and 
Baltimore, Maryland, metropolitan area is examined by incorporating trees, soil, and grass into the coupled 
Weather Research and Forecasting model and an urban canopy model (WRF-UCM). By parameterizing the 
effects of these natural surfaces alongside roadways and buildings, the modified WRF-UCM is used to in- 
vestigate how urban trees, soil, and grass dampen the UHI. The modified model was run with 50% tree cover 
over urban roads and a 10% decrease in the width of urban streets to make space for soil and grass alongside 
the roads and buildings. Results show that, averaged over all urban areas, the added vegetation decreases 
surface air temperature in urban street canyons by 4.1 K and road-surface and building-wall temperatures by 
15.4 and 8.9 K, respectively, as a result of tree shading and evapotranspiration. These temperature changes 
propagate downwind and alter the temperature gradient associated with the Chesapeake Bay breeze and, 
therefore, alter the strength of the bay breeze. The impact of building height on the UHI shows that decreasing 
commercial building heights by 8 m and residential building heights by 2.5 m results in up to 0.4-K higher 
daytime surface and near-surface air temperatures because of less building shading and up to 1.2-K lower 
nighttime temperatures because of less longwave radiative trapping in urban street canyons. 


1. Introduction 

Urbanization can alter local climate and form an ur- 
ban heat island (UHI; Landsberg 1981). Altering land 
use by creating impervious urban surfaces causes in- 
creased runoff, decreased evapotranspiration, increased 
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solar radiation absorption, additional release of an- 
thropogenic heat, and changes in surface friction, which 
results in changes in near-surface air temperature, hu- 
midity, wind speeds, low-level convergence/divergence, 
convection, and precipitation (e.g., Oke and Cleugh 
1987; Bornstein and Lin 2000; Arnfield 2003). Previous 
studies show that UHIs strengthen as city size and 
building density increase (Oke 1973; Landsberg 1981; 
Atkinson 2003; Imhoff et al. 2010). In addition, a recent 
study reveals that upstream urbanization can magnify 
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UHI effects under the favorable influence of larger-scale 
flow (Zhang et al. 2009, 2011). Therefore, a smaller city 
with upstream urbanization can have a larger UHI than 
a larger city with no upstream urbanization (Zhang et al. 
2009). 

A UHI can have ill effects on human health. The UHI 
can amplify summertime heat waves, leading to heat 
stress (Kunkel et al. 1996). The UHI can also aggravate 
air pollution. Air quality model results show air pollu- 
tion worsens as temperature increases (Weaver et al. 
2009; Banta et al. 1998; Cheng and Byun 2008; Jacob and 
Winner 2009), and observations confirm a correlation 
between high temperature and events with high air 
pollution (Bloomer et al. 2009, 2010; Tai et al. 2010). 

Urban trees have the potential to dampen the UHI 
and decrease near-surface temperatures through direct 
shading and evaporative cooling. Observations show 
temperatures over a grassy surface were 0.7-1. 3 K 
cooler under urban trees than adjacent areas with no 
tree cover (Souch and Souch 1993). Similar measure- 
ments revealed temperatures to be 2. 2-3. 3 K cooler 
under mature trees in a suburban neighborhood than 
in new developments with no trees (McGinn 1982). 
Planting trees near buildings can reduce summertime 
building-surface temperatures through shading and 
evapotranspiration, which can in turn reduce energy 
use for cooling building interiors (Heisler 1986) and 
therefore emissions from power plants. Specifically, 
trees planted near a building reduce building heat 
gain by 1) shading solar radiation incident on the 
building, 2) shading nearby surfaces that radiate heat 
toward the building, 3) reducing the outside air in- 
filtration rate by lowering ambient wind speeds, and 4) 
lowering ambient temperatures through evapotrans- 
piration (Akbari 2002). A summertime energy-use 
study in which eight 6-m-tall and eight 2.4-m-tall trees 
were planted around one house and then moved to 
another house showed that shading by these trees de- 
creased summertime energy demand by as much as 
30% (Akbari et al. 1997). 

Planting trees in urban areas can also improve air 
quality by decreasing air temperatures and increasing 
the removal rate of pollutants from the atmosphere. 
A decrease in air temperature lowers the temperature- 
dependent emission rate of volatile organic compounds 
(VOCs) from vegetation and temperature-dependent 
chemical reaction rates that produce tropospheric ozone 
(Taha 1995). Although increasing the amount of trees 
increases the amount of vegetation that emits VOCs 
that react with nitrogen oxides to form ozone, trees also 
filter air pollution. Most gaseous air pollutants that are 
removed by trees are taken up by leaf stomata, but some 
are removed by the surface of trees (Nowak et al. 2006). 


Some particles that come into contact with trees are 
absorbed into the trees, but most remain on the surface 
of trees until they are resuspended into the atmosphere 
or are transported to the surface by precipitation or 
falling leaves or branches (Nowak et al. 2006). A re- 
gional meteorological and air quality modeling study 
revealed that planting trees that emit low amounts of 
VOCs decreased tropospheric ozone concentrations 
whereas planting trees that emit large amounts of VOCs 
degraded air quality in Los Angeles, California (Taha 
1995). This study did not account for further reductions in 
ozone concentrations due to changes in energy demand 
and anthropogenic emissions from planting additional 
trees. A further modeling study shows that planting 20 
million trees in the Los Angeles air basin reduces the 
mass of ozone in the mixed layer by 4.5% because of 
increased deposition of ozone and nitrogen dioxide (Taha 
et al. 1997). 

A UHI in a coastal region may have the potential to 
alter a sea, lake, or bay breeze, thereby affecting air 
quality. A summertime sea breeze is driven by the tem- 
perature gradient between the warm land and cool water. 
Increasing or decreasing land surface temperature through 
urbanization or planting urban trees, respectively, will al- 
ter the temperature gradient and therefore modify the 
strength of the sea breeze. Previous studies have shown 
that a sea-breeze circulation can exacerbate air pollution 
levels (Boucouvala and Bornstein 2003; Evtyugina et al. 
2006; Loughner et al. 2011). In Houston, high ozone epi- 
sodes begin when the large-scale flow is offshore before 
a bay breeze develops (Banta et al. 2005; Darby 2005). As 
the bay breeze begins to strengthen, stagnant conditions 
develop, allowing ozone and its precursors to accumulate 
before being advected farther onshore as the bay breeze 
increases in intensity later in the afternoon (Banta et al. 
2005; Darby 2005). 

With the rapid increases in computing power in recent 
years, there have been a growing number of higher- 
resolution model simulations. Many studies show bene- 
fits of using high-resolution mesoscale models to resolve 
frontal structures, orographical flows, and vertical circu- 
lations induced by surface inhomogeneities [see Mass 
et al. (2002) for a review]. For urban settings at finescales 
(horizontal grid spacing less than 1 km), urban canopy 
models are used to simulate the meteorological conditions 
in the complex urban environment consisting of streets, 
buildings, and vegetation. Currently, the Weather Re- 
search and Forecasting (WRF) model (Skamarock et al. 
2008) can be run coupled with the “Noah” land surface 
model and an urban canopy model (Kusaka et al. 2001) 
(WRF-UCM). This urban canopy model does not include 
soil, grass, or trees in urban street canyons, however. Lee 
and Park (2008) developed a vegetated urban canopy 
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model and found that canyon vegetation has a large in- 
fluence on surface temperatures and sensible and latent 
heat fluxes. 

In the present study, the effects of soil, grass, and trees 
in urban street canyons are incorporated into version 
3.1.1 of the WRF-UCM to examine how urban trees 
dampen the UHI and how sensitive the UHI is to urban 
building height. These objectives are accomplished by 
performing simulations down to a horizontal grid spac- 
ing of 0.5 km with a modified version of the Advanced 
Research WRF-UCM in which urban trees, soil, and 
grass are incorporated, following closely the work of Lee 
and Park (2008). The simulations cover the Washington- 
Baltimore metropolitan areas from 1200 UTC 7 July 
to 1200 UTC 10 July 2007. High air pollution was ob- 
served in the region under hot, sunny, stagnant con- 
ditions. A cold front passed through the area on 6 July 
2007, followed by a short-wave trough approaching the 
Mid- Atlantic region (Zhang et al. 2011). The Baltimore, 
Maryland, UHI exhibited a 2-m temperature of 37.5°C, 
and downwind of Baltimore near-surface 8-h maximum 
ozone mixing ratios reached 125 ppb (the current air 
quality standard is 75 ppb) on 9 July. Previous work ex- 
amining this modeling scenario explored the impact of 
upstream urbanization on the UHI (Zhang et al. 2009, 
2011), characterized the air pollution event (Yegorova 
et al. 2011), and investigated the effects of fair-weather 
cumulus clouds and the Chesapeake Bay breeze on air 
quality (Loughner et al. 2011). 

2. Model description and modification 

In this section, we will describe the basic configura- 
tion of the coupled WRF-UCM model and then show 
how the tree effects are incorporated into the coupled 
model, followed by the description of sensitivity- 
experiment designs to study the impact of trees on 
UHI. 

a. Model configuration 

In the coupled Advanced Research WRF-UCM, the 
Noah land surface model calculates soil moisture and 
temperature, skin temperature, canopy water content, 
and the energy and water flux terms in the surface and 
water energy balance equations (Chen and Dudhia 2001) 
and the UCM improves the parameterization of physical 
processes involved in the exchange of heat, momentum, 
and water vapor in urban environments by including 
shadowing from buildings, reflection of shortwave and 
longwave radiation, wind profile information in the 
canopy layer, and a multilayer heat transfer equation 
for roof, wall, and road surfaces (Kusaka and Kimura 
2004). 


The WRF-UCM includes three categories of ur- 
ban surfaces: commercial/industrial/transportation, 

high-intensity residential, and low-intensity residential. 
Each urban category consists of fractional coverage of 
urban land occupied by buildings and roads, with the 
remaining fraction as undeveloped land. Urban fraction is 
set to 95%, 90%, and 50% for grid cells labeled com- 
mercial/industrial/transportation, high-intensity residen- 
tial, and low-intensity residential, respectively, for the 
model simulations described herein. The remaining land 
in grid cells classified as urban is considered to be un- 
developed and is classified as U.S. Geological Survey 
(USGS) land-use type “cropland/grassland mosaic.” Sur- 
face heat and moisture fluxes are calculated in the UCM 
for urban streets and buildings and in the Noah land sur- 
face model for undeveloped urban land. 

The coupled model is run at 13. 5-, 4.5-, 1.5-, and 
0.5-km horizontal grid spacings with x, y dimensions of 
181 X 151, 244 X 196, 280 X 247, and 349 X 349 grid 
cells, respectively (see Fig. 1 for the model domains). All 
of the domains use 32 layers in the vertical direction with 
20 layers in the lowest 2 km. The North American 
Regional Reanalysis is used for the model initial and 
outermost lateral boundary conditions. The three cat- 
egories of urban surfaces are defined by the U.S. En- 
vironmental Protection Agency’s National Land Cover 
Dataset for 2001. The urban areas in the 0.5-km domain 
are displayed in Fig. 2. The model physics schemes used 
include 1) a double-moment six-class microphysics 
scheme (Lim and Hong 2010), 2) the Mellor-Yamada- 
Janjic boundary layer parameterization (Janjic 1994), 
3) the Noah land surface model (Chen and Dudhia 

2001) , and 4) an ensemble cumulus parameterization 
that advances the Grell-Devenyi scheme (Grell and 
Devenyi 2002) to allow subsidence in neighboring 
grid cells (Skamarock et al. 2008). The cumulus pa- 
rameterization is only used for the 4.5- and 13.5-km 
domains. 

b. Including trees in the UCM 

The WRF-UCM was modified to investigate how ur- 
ban trees impact the UHI effects. The WRF-UCM was 
tailored to include grass, soil, and trees in urban street 
canyons and trees in undeveloped land in urban grid 
cells. Urban vegetation coverage is described in Table 1. 
The urban tree leaf area index (LAI) is set to 3.31, the 
same as USGS land-use type “deciduous broadleaf 
forest” and in line with favorable urban trees found in 
the Mid-Atlantic region. Averaged LAIs by individual 
tree species from a tree survey in New York City, New 
York, range from 0.68 for honey locust trees to 10.07 for 
northern white cedar trees (TreesNY and CENYC 

2002) . Ideally, urban trees must be durable, must be able 
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Fig. 1. Location of model domains 1, 2, 3, and 4, which have horizontal grid spacings of 13.5, 4.5, 

1.5, and 0.5 km, respectively. 


to thrive in an urban environment, must be low biogenic 
VOC emitters, and must have a large LAI to cause 
maximum shading, latent heat exchange, and pollutant 
deposition. A list of urban tree species from the U.S. 
Department of Agriculture Forest Service (Chicago 
Botanic Garden 2011) was referenced to find trees that 
are suitable to live in an urban environment. Biogenic 
emissions of tree species from the Biogenic Emissions 
Inventory System (Vukovich and Pierce 2002) were 
referenced to find low-VOC emitters. LAIs from a New 
York City tree survey (TreesNY and CENYC 2002) 
were referenced to determine which trees have a large 
LAI. It was determined that ash and elm trees are pre- 
ferred urban trees for the Washington and Baltimore 
metropolitan areas. Observed average LAI in New 
York City for ash and elm trees is 4.11 and 3.12, re- 
spectively (TreesNY and CENYC 2002). Urban-canyon 
tree height is set to 10 m, which is in line with observed 
ash and elm tree heights in New York City. Average 
and maximum tree heights in New York City are 11.84 
and 19.50 m, respectively, for ash trees and 11.58 and 
27.40 m for elm trees (TreesNY and CENYC 2002). 

To increase the percentage tree cover in the un- 
developed land from 0% to 50%, the undeveloped land 
is reclassified from USGS land-use type cropland/ 
grassland mosaic to USGS land-use type “cropland/ 
woodland mosaic.’’ A 50% tree cover over undeveloped 
land and urban streets results in 26.25%, 27.5%, and 37.5% 
tree cover over commercial/industrial/transportation, high- 
intensity residential, and low-intensity residential urban 
areas, respectively. Tree canopy cover in urban and met- 
ropolitan areas in the United States averages 27% and 


33%, respectively (Dwyer and Nowak 2000). Surface heat 
and moisture fluxes are calculated with the Noah land 
surface model for undeveloped urban land and with the 
UCM for the urban buildings and street canyons. 

The UCM was modified to account for fractional 
coverage of grass, soil, and tree canopies in the street 
canyons, including the increased momentum drag due to 
the tree canopy; the transmissivity of shortwave and 
longwave radiation through the tree canopy; tree shad- 
ing on building roofs, building walls, and the ground; and 
additional shortwave radiative, longwave radiative, la- 
tent heat, and sensible heat fluxes due to the added trees, 
soil, and grass. The UCM was modified by blending 
a vegetated urban canopy model (Lee and Park 2008) 
with a single-layer urban canopy model (Kusaka et al. 
2001; Kusaka and Kimura 2004), which is already cou- 
pled inside the WRF modeling system. Details on the 
blending of these two models are described in the ap- 
pendix. 

c. Experimental design 

Three different WRF-UCM simulations were per- 
formed to achieve the above-mentioned objectives: 1) 
a base case; 2) a run with urban soil, grass, and trees; and 
3) a scenario with shorter buildings. Hereinafter, these 
simulations are referred to as “no trees,” “trees,” and 
“no trees/shorter buildings.” Urban fraction and build- 
ing and canyon dimensions for the three model simu- 
lations and three urban categories are displayed in 
Table 2. The no trees/shorter buildings simulation is per- 
formed to test the sensitivity of the UHI effects to building 
height. 
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Fig. 2. Urban land use in the 0.5-km-horizontal-resolution domain, where red, yellow, and 
cyan represent commercial/industrial/transportation, high-intensity residential, and low- 
intensity residential, respectively. The letters F, R, B. and W denote the cities of Frederick, MD; 
Reston, VA; Baltimore; and Washington. The letter C shows the location of the Chesapeake 
Bay. 


3. Sensitivity simulations 

Three nested-grid simulations (i.e., no trees, trees, and 
no trees/shorter buildings) with the finest grid spacing of 
0.5 km are analyzed alongside observations to in- 
vestigate the role of urban trees and building heights on 
the UHI. Averaged over 23 measurement sites within 
the 0.5-km-resolution domain, the near-surface (at 
height z = 2 m) temperatures for all three simulations 
are compared with temperature and wind velocity ob- 
servations from the National Weather Service and the 
Maryland Department of the Environment. The obser- 
vational sites include urban sites and sites that are 
downwind and upwind of urban areas. The model sites 


include the observational sites but are averages over 
a 0.5-km domain that may include a mixture of urban 
streets, buildings, and undeveloped land. Statistics of 
2-m temperature comparisons between the model sim- 
ulations and observations show that all three model 
simulations agree reasonably well with the observations 
(Table 3). There are noticeable differences among the 
simulations, however. The trees simulation consistently 
has lower surface temperatures than the no-trees simu- 
lation because of increased evapotranspiration and 
shading of roads and buildings. At the measurement 
sites, the surface temperature of the trees simulation is 
on average 0.3 K cooler than the no-trees simulation 
between 1200 UTC 8 July and 1200 UTC 10 July. The no 


Table 1. Specification of urban-canyon grass/soil fraction, urban-canyon tree canopy fraction, undeveloped-land tree canopy fraction, 
and tree cover fraction over total grid cells for the three urban categories [low-intensity residential (LI), high-intensity residential (HI), 
and commercial/industrial/transportation (C)] and the three simulations [no trees (NT), trees (T), and no trees/shorter buildings (NTSB)]. 



Urban-canyon grass/ 
soil fraction 

Urban-canyon tree 
canopy fraction 

Undeveloped-land 
tree canopy fraction 

Tree cover of 
total grid (%) 

LI 

HI 

c 

LI 

HI 

c 

LI 

HI 

c 

LI 

HI 

c 

NT 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

T 

0.1 

0.1 

0.1 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

37.5 

27.5 

26.25 

NTSB 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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TABLE 2. Specification of urban fraction, building height and depth, and urban-canyon width for the three urban categories and three 

simulations. 



Urban fraction 

Building height (m) 

Building depth (m) 

Urban-canyon width (m) 

LI 

HI 

C 

LI 

HI 

C 

LI 

HI 

C 

LI 

HI 

C 

NT 

0.5 

0.9 

0.95 

7.5 

10 

20 

8.3 

9.4 

10 

8.3 

9.4 

10 

T 

0.5 

0.9 

0.95 

7.5 

10 

20 

8.3 

9.4 

10 

8.3 

9.4 

10 

NTSB 

0.5 

0.9 

0.95 

5 

7.5 

12 

8.3 

9.4 

10 

8.3 

9.4 

10 


trees/shorter buildings simulation has higher daytime 
temperatures but lower nighttime temperatures than the 
no-trees simulation. Averaged over all of the measure- 
ment sites, maximum differences between the no trees/ 
shorter buildings and no-trees simulations peaked at 0.4 K 
during the day and 1.2 K during the night. Shorter build- 
ings produce fewer shadows, thereby allowing more solar 
radiation to heat the building walls and roads. On the 
other hand, street canyons with shorter buildings trap less 
longwave radiation emitted from the surfaces, allowing 
the surface to cool more quickly during the nighttime. 

At 2000 UTC (i.e., 1500 LST) 9 July 2007, the Ches- 
apeake Bay breeze’s convergence zone was over Balti- 
more (Fig. 3), and at 2300 UTC (i.e., 1800 LST) it was 
between Washington and Baltimore (Fig. 4). The no 
trees/shorter buildings simulation is slightly warmer than 
the no-trees simulation. The maximum 2-m temperature 
difference between the no trees/shorter buildings and no- 
trees simulations in Washington is 0.6 K at 2000 UTC and 
0.4 K at 2300 UTC 9 July (not shown). Even though this 
difference is small, it can be seen from Figs. 3 and 4 that 
the spatial extent of the highest temperatures in the re- 
gion is largest in the no trees/shorter buildings simulation, 
followed by the no-trees and then the trees simulations. 

The maximum difference between the no-trees- and 
trees-simulated 2-m temperature is 0.8 and 1.2 K at 2000 
and 2300 UTC, respectively, in Washington and 1.9 K at 
2000 and 2300 UTC in Baltimore. Differences in 2-m 
temperature between the no-trees and the trees simu- 
lations at 2000 (Fig. 5a) and 2300 UTC (Fig. 5b) show 
that the impact of planting trees in all urban areas in the 
domain is largest in Baltimore. 

Urban trees are seen to affect the strength of the 
Chesapeake Bay breeze. Figures 5a and 5b show the 
change in temperature associated with trees at 2000 and 
2300 UTC, respectively. The most striking feature is a 
thin line of warmer temperatures in the trees simulation 
than in the no-trees simulation running north-south 
between Washington and Baltimore. Northwest of Bal- 
timore, this thin line switches sign, indicating cooler 
temperatures in the trees simulation. At 2300 UTC, the 
thin line shows temperatures of up to 2 K warmer in the 
trees simulation than in the no-trees simulation east and 
northeast of Washington but temperatures that are 2 K 


cooler northwest of Baltimore (Fig. 5b). This thin line is 
attributable to differences in the positioning of the 
Chesapeake Bay-breeze convergence zone. Fligh tem- 
peratures resulting from the Washington UFII propagat- 
ing downwind near the coastline of the Chesapeake Bay 
strengthen the temperature gradient along the Ches- 
apeake Bay coastline and therefore strengthen the bay 
breeze. The addition of urban trees dampens the urban 
heat island, causing a weaker temperature gradient along 
the coastline and a weaker bay breeze that does not 
penetrate as far inland downwind of Washington. So, the 
thin line of warmer temperatures in the trees simulation 
shown in Figs. 5a, 5b, and 5d to the east and northeast of 
Washington can be attributed to a weaker bay breeze that 
does not push the cooler air originating over the water as 
far inland. On the other hand, once a bay breeze pene- 
trates inland past an urban area with cool air originating 
over the water moving over the city, the UFII warms this 
cool air, which is being transported to the bay-breeze 
front. This warming of the cool air causes a weaker tem- 
perature gradient along the bay-breeze front, lessening 
the inland penetration of the bay breeze. The addition of 
trees to a coastal city, however, causes dampening of the 
UHI so that the amount of warming of the cool air moving 
over the city toward the inland bay-breeze front by the 
UHI is reduced. Therefore, for conditions in which a bay 
breeze penetrates past a coastal city, the addition of trees 
results in a stronger bay breeze that penetrates farther 
inland. So, the thin line of cooler temperatures in the trees 
simulation shown in Figs. 5a, 5b, and 5c to the northwest 
of Baltimore can be attributed to a stronger bay breeze 
that penetrates farther inland. As mentioned previously, 


TABLE 3. Mean absolute error (MAE), mean bias error (MBE), 
and root-mean-square error (RMSE) of 2-m temperature from 
0000 UTC 8 Jul through 1200 UTC 10 Jul 2007 for the three sim- 
ulations calculated at 23 measurement sites within the 0.5-km 
horizontal domain. The observational locations include sites lo- 
cated in, downwind, and upwind of urban areas. 



MAE 

MBE 

RMSE 

NT 

1.75 

0.50 

2.23 

T 

1.74 

0.24 

2.25 

NTSB 

1.65 

0.41 

2.10 
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a) Observations 



b) No Trees Shorter Buildinqs 
39. 6N 


39. IN 



38. 6N 

77.35W 76.85W 76.35W 77.35W 76.85W 76.35W 


c) Trees d) No Trees 



2 m temperature (K) 10 m wind speed 


Fig. 3. Observed and simulated 2-m temperature (color shading) and 10-m wind speed (arrows) at a horizontal 
resolution of 0.5 km at 2000 UTC (i.e., 1500 LST) 9 Jul 2007. 


the impact of planting trees in all urban areas in the do- 
main is larger in Baltimore than in Washington. Since the 
bay breeze is stronger over and downwind of Baltimore in 
the trees simulation than in the no-trees simulation after 
the bay-breeze convergence zone passes through Balti- 
more, more air originating over the cool surface waters is 
transported to Baltimore. So, cooler temperatures over 
Baltimore in the trees simulation than in the no-trees 
simulation are attributable to a stronger bay breeze 
transporting cool air over Baltimore at a faster rate. 


There are significant urban temperature differences 
on the subgrid scale. Figure 6 shows building-roof-, 
building-wall-, and road-surface temperatures as well as 
street-canyon air temperature averaged over all urban 
buildings, roads, and street canyons for the three simu- 
lations. The roof-surface temperatures vary little among 
the three simulations because the roofs for both the no- 
trees and no trees/shorter buildings simulations are not 
shaded and only the low-intensity-residential roofs in 
the trees simulations are partially shaded at low solar 
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a) Observations 


b) No Trees Shorter Buildings 



c) Trees 
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FIG. 4. As in Fig. 3, but at 2300 UTC (i.e., 1800 LST). 


zenith angles. The no trees/shorter buildings simulation’s 
maximum daytime building-wall- and road-surface tem- 
peratures and maximum canyon air temperature are 1.9, 
2.4, and 1.5 K greater, respectively, than the no-trees 
simulation’s temperatures on 9 July because of more 
shading from the buildings. On the other hand, the no 
trees/shorter buildings simulation’s minimum nighttime 
building-wall- and road-surface temperatures and mini- 
mum canyon air temperature are 0.4, 0.6, and 0.6 K 
lower, respectively, than in the no-trees simulation’s 
temperatures between sunset on 8 July and sunrise on 


9 July as a result of less longwave radiative trapping in 
the urban canyon. The trees simulation’s maximum 
building-wall- and road-surface temperatures and max- 
imum canyon air temperature are 8.9, 15.4, and 4.1 K 
lower, respectively, than the no-trees simulation’s 
temperatures on 9 July because of tree shading and 
evapotranspiration. The trees simulation’s minimum 
wall- and road-surface temperatures and canyon air tem- 
perature are 3.1, 3.2, and 2.5 K lower, respectively, than 
the no-trees simulation’s temperatures during the night 
of 8 July. 
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a) 2000 UTC 9 July 2007 b) 2300 UTC 9 July 2007 



2007 


39. 6N 


39.55N- 



39. 5N 


- 1.5 -1 - 0.5 0.5 1 1.5 2 

2 m temperature difference (K) 


76. 8W ? ^-,76.75W 
10 m wind speed 
difference 


Fig. 5. Difference in trees minus no-trees 2-m temperature simulations with a horizontal resolution of 0.5 km at (a) 
2000 UTC (i.e., 1500 LST) and (b) 2300 UTC (i.e., 1800 LST) 9 Jul 2007. The areas in the rectangles in (b) that are 
northwest of Baltimore and between Washington and Baltimore are shown in detail in (c) and (d), respectively, for 
2300 UTC, with 10-m wind velocities added. 


4. Concluding remarks 

In this study, the roles of urban trees and building 
height in UHI effects are investigated by performing 
two simulations with a modified version of the WRF- 
UCM: one with the inclusion of the parameterized ef- 
fects of urban trees, soil, and grass and the other without 
those effects. Results show that urban areas with shorter 
buildings have a larger diurnal cycle. Shorter urban 


buildings cause higher surface and near-surface air 
temperatures during the daytime because of less build- 
ing shading and cause lower temperatures at night as 
a result of less longwave radiative trapping in urban 
street canyons. In the model simulations, decreasing 
building size causes maximum daytime urban-canyon air 
temperature to increase by 1.5 K and minimum night- 
time urban-canyon air temperature to decrease by 0.6 K. 
Urban trees result in lower surface and near-surface 
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a) Roof Surface 



12UTC 7 July 12UTC 8 July 12UTC 9 July 
Time and Date 


b) Building Wall Surface 



12UTC 7 July 12UTC 8 July 12UTC 9 July 1 2UTC 10 Ju 
Time and Date 



1 2UTC 7 July 1 2UTC 8 July 1 2UTC 9 July 
Time and Date 


d) Canyon Air 



12UTC 7 July 12UTC 8 July 12UTC 9 July 12UTC 10 Ju 
Time and Date 


Fig. 6. Time series of no trees/shorter buildings minus no trees (blue), and trees minus no trees (green) simulated 
subgrid (a) roof-, (b) building-wall-, and (c) road-surface temperatures and (d) canyon air temperature with a hori- 
zontal resolution of 0.5 km averaged over all urban land-use categories from the 3-day simulations between 
1200 UTC 7 Jul and 1200 UTC 10 Jul 2007. 


air temperatures because of tree shading and evapo- 
transpiration. The addition of trees in the simulations 
causes maximum daytime and minimum nighttime urban- 
canyon air temperatures to decrease by 4.1 and 2.5 K, 
respectively. Future investigations will look into how 
these temperature differences that result from varying 
building height and urban vegetation cover influence 
emissions, climate, and air quality. 

Urban trees alter the strength of the Chesapeake Bay 
breeze by altering the temperature gradient near the 
coastline. Urban trees on the warm side of the Ches- 
apeake Bay breeze dampen the strength of the breeze 
and therefore suppress the penetration of the bay breeze 
inland. The urban trees decrease the near-surface air 
temperature over the warm land and therefore decrease 
the temperature gradient between the warm air over 
land and relatively cool air over the water. For a sce- 
nario in which the bay breeze penetrates inland past 
Baltimore, however, it is found that additional urban 
trees in Baltimore increase the strength of the bay 
breeze. The addition of urban trees in Baltimore causes 
the temperature of the cool side of the bay breeze to 
decrease, resulting in a stronger temperature gradient 


and bay breeze. Because the bay breeze can have an 
impact on air quality, future research can examine how 
altering the bay breeze through the addition of urban 
trees can affect air quality. 

Because urban areas in the WRF-UCM only include 
impervious surfaces, buildings, and roads, the WRF- 
UCM with trees can be used to simulate the meteoro- 
logical conditions under more realistic urban land surface 
coverage. The National Land Cover dataset used in this 
study is from 2001. Keeping urban tree fraction and 
building height datasets up to date is essential for ini- 
tializing the WRF-UCM with trees, given that the UHI 
is sensitive to these parameters. 

WRF-UCM with trees can be used as a tool to in- 
vestigate how vegetative land surface changes to urban 
areas affect many aspects of earth science. Results show 
that urban trees can lower surface and air temperatures 
in and downwind of cities and alter a bay breeze. Future 
studies will investigate how lower summertime building- 
surface temperatures due to an increase in urban trees 
can result in less energy demand for cooling buildings, 
will quantify the reduction in energy demand to a de- 
crease in anthropogenic emissions, and will determine 


October 2012 


LOUGHNER ET AL. 


1785 


the impact of urban trees on climate and air quality. The 
WRF-UCM with trees can be utilized to investigate how 
changes in urban trees, grass, and/or soil affect the hy- 
drological cycle through changes in runoff, evapotrans- 
piration, and even precipitation both in and downwind of 
cities. 
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APPENDIX 


LAl a = (h f - h b )p,, (Al) 

and the LAI in the street canyon is defined as 

LAI c = Vg ( A2 ) 

where hf and h b are the tree and building heights, re- 
spectively. The density of the leaves is assumed to be 
invariant with height and is defined by 

Pl = LAllh f . (A3) 

Table Al lists all of the variables used to parameterize 
urban trees in the WRF-UCM. 

The amount of radiation that can be transmitted 
through a tree canopy needs to be defined to compute 
the radiative budget. The transmissivity of radiation 
through a tree canopy from height z\ to z 2 is computed by 


Parameterization of Trees in the UCM 

Energy storage and surface temperatures are calcu- 
lated on the basis of the surface energy balance equation 
described in both Lee and Park (2008) and Kusaka et al. 
(2001). The UCM is coupled with the Noah land surface 
model to obtain the latent and sensible heat fluxes at the 
soil surface by passing the incident longwave and 
shortwave radiation at the soil surface. The sensible heat 
flux of the tree canopy and evapotranspiration flux, 
which include the moisture flux on tree leaves and 
transpiration from the root zone, are calculated follow- 
ing Lee and Park (2008). Transpiration from the root 
zone is a function of stomatal resistance, which is also 
obtained from the Noah land surface model. Wind speed 
below treetops is modified to account for the additional 
trees by using the method of Lee and Park (2008). 

The longwave and shortwave radiation fluxes are 
calculated by blending a vegetated urban canopy model 
(Lee and Park 2008) and an equation for transmission of 
radiation through a tree canopy (Annandale et al. 2004; 
Campbell and Norman 1998; Norman and Welles 1983) 
with the single-layer urban canopy model (Kusaka et al. 
2001; Kusaka and Kimura 2004), which is coupled within 
the WRF modeling system. For completeness, we review 
many of the equations used in the vegetated urban 
canopy model and the single-layer urban canopy model 
before we present the final blended longwave and 
shortwave radiation flux equations on the surfaces of 
buildings, roads, soil/grass, and trees to introduce the 
variables used in the final merged equations. 

For trees taller than the adjacent buildings, the tree 
canopy is split into two layers: one above the building 
roof and another in the street canyon. The LAI above 
the street canyon is calculated by 


T(z v z 2 ) = exp 


—0.5 p l \fa dz 


(A4) 


where p/ is the leaf area density and a represents the 
absorptivity for individual leaves (Annandale et al. 2004; 
Campbell and Norman 1998; Norman and Welles 1983). 
A typical value of a is 0.5 (Annandale et al. 2004) and is 
used here. 

The direct solar radiation that reaches the surface of 
the leaves in the street canyon is 


Sr=S? l f l {\\-T(h c ,h {c )]sme n 

Ijl-mya-siny}, (A5) 

where Sf* h c , and hf C are the downward direct solar 
radiation that reaches the surface of leaves, downward 
direct solar radiation in the canyon, the height of the top 
of the tree in the canyon, and the average tree shaded 
height due to the buildings, respectively (Lee and Park 
2008); here, 0 n is the angle between the average canyon 
axis and the direction of the sun as described by Kusaka 
et al. (2001). The first and second terms on the rhs of 
Eq. (A5) describe the solar radiation that reaches 
building shaded and unshaded leaves, respectively. If 
the tree height is greater than the adjacent building 
height, then h [c is set to the building height. The height 
of the shaded canopy is defined by 


h = h 


fc 


W 

2 tan 0 Z sin Q n ’ 


(A6) 


where w is the width of the ground between the buildings 
and 0 Z is the solar zenith angle. If h c is greater than h tc , 
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Table Al. Description of symbols used to parameterize trees in UCM. 

Symbol 

Description 

Units 

LAI 

Leaf area index 

— 


LAI* 

Effective leaf area index 

— 


T, 

Temperature of leaves 

K 


Tl a 

Temperature of leaves above roof level 

K 


Y 

Extinction coefficient 

— 


£/ 

Emissivity of leaves 

— 


pi 

Leaf density 

nU 1 


(j 

Leaf aspect ratio 

— 


a 

Absorptivity of a leaf 

— 


f, 

Tree canopy fraction in street canyon 

— 


/rd 

Road fraction of ground in street canyon 

— 


fs 

Grass/soil fraction of ground in street canyon 

— 


h f 

Tree height 

m 


hb 

Building height 

m 


h c 

Height of top of tree in street canyon 

m 


hfc 

Average tree-shaded height due to buildings 

m 


hf a 

Height between the treetop and building roof level 

m 


ha—>r direct 

Distance of path of direct radiation through the tree canopy from the center of the shadow on roof 

m 


e„ 

Angle between average street axis and sun angle 

rad 


e z 

Solar zenith angle 

rad 


a r 

Albedo of roof 

— 


w 

Width of the ground between the buildings 

m 


^shadow 

Normalized shadow length 

— 


^height 

Normalized building height 

— 


Aground 

Normalized street-canyon width 

— 


^shadow roof 

Normalized shadow on roof 

— 


Hzu zi) 

Transmissivity between heights Zi and z 2 

— 


T wa 

Wall-atmosphere transmissivity 

— 


T ww 

Wall-wall transmissivity 

— 


T wg 

Wall-ground transmissivity 

— 


T ga 

Ground-atmosphere transmissivity 

— 


S D* 

Downward direct solar radiation reaching leaves 

W nY 


s : f 

Downward direct solar radiation in street canyon 

W nT 

-2 

cD l 

Downward direct solar radiation reaching wall 

W nY 

-2 

oD* 

^rd 

Downward direct solar radiation reaching road 

W nY 

-2 

s f l 

Downward direct solar radiation reaching grass/soil 

W nY 

-2 

s;‘ 

Indirect solar radiation absorbed by leaves 

W nY 

-2 

s{ 1 

Indirect solar radiation in street canyon 

W nY 

-2 

s 1 * 

•Y 

Indirect solar radiation reaching leaves if no buildings present 

W nY 

-2 

S C 

Indirect solar radiation incident on wall 

W nY 

-2 

S 1 J 
J rd 

Indirect solar radiation incident on road 

W nY 

-2 

5 71 

Indirect solar radiation incident on grass/soil 

W nY 

-2 

511 

Solar radiation absorbed by wall 

W nY 

-2 

^d 1 

Solar radiation absorbed by road 

W nY 

-2 

s s n 

Solar radiation absorbed by grass/soil 

W nY 

-2 

ct i 

Solar radiation absorbed by leaves 

W nY 

-2 

s o* 

Direct solar radiation that reaches leaves above roof level 

W nY 

-2 

s?‘ 

Direct solar radiation from the atmosphere above 

W nY 

-2 

< 

Indirect solar radiation that reaches leaves above roof level 

W nY 

-2 

Si 1 

Indirect solar radiation from the atmosphere above 

W nY 

-2 

oD* 
15 la 

Direct solar radiation absorbed by leaves above roof level 

W nY 

-2 

o/i 

"’la 

Indirect solar radiation absorbed by leaves above roof level 

W nY 

-2 

oD^ 

Direct solar radiation absorbed through side of tree canopy above roof level 

W nY 

-2 

o/l 

J la_side 

Indirect solar radiation absorbed through side of tree canopy above roof level 

W nY 

-2 

Solar radiation reflected by roof and absorbed by tree canopy above roof 

W nY 

-2 

5U 

Solar radiation absorbed by roof 

W nY 

-2 

L cl 

Longwave radiation entering street canyon from aloft 

W nY 

-2 

L\ t 

Longwave radiation emitted by leaves that reaches other leaves 

W nY 
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TABLE Al. (Continued) 


Symbol Description Units 


Longwave radiation emitted by leaves that escapes canyon W m~ 2 

LJiJ Longwave radiation emitted by leaves that reaches wall W m~ 2 

Longwave radiation emitted by leaves that reaches ground W m~ 2 

Lll Net longwave radiation on surface of wall W m~ 2 

L'J Net longwave radiation on surface of road W m~ 2 

Lj 1 Net longwave radiation on surface of grass/soil W m~ 2 

L/* Net longwave radiation on surface of leaves W m~ 2 

L l l2 Longwave radiation emitted by leaves that is absorbed by leaves W m~ 2 

L( d2 Longwave radiation emitted by the road that is absorbed by leaves W m~ 2 

L \ 2 Longwave radiation emitted by soil that is absorbed by leaves W m~ 2 

L l w2 Longwave radiation emitted by building walls that is absorbed by leaves W m~ 2 

L[ 2 Downwelling longwave radiation from above the urban canopy that is absorbed by leaves W m~ 2 

Longwave radiation emitted from the tree canopy above the rooftop that enters the street canyon W m~ 2 

Li aT Longwave radiation emitted from leaves above roof level W m~ 2 

L|“ Longwave radiation emitted from leaves that reaches the surface of other leaves above the roof height W m~ 2 

Ly Longwave radiation absorbed by roof W m~ 2 

Longwave radiation absorbed by leaves above roof level W m~ 2 


F w ^a Wall-atmosphere view factor 

F g „ a Ground-atmosphere view factor 

F w ^ s Wall-ground view factor 

F w „„ Wall-wall view factor 

Leaves-atmosphere view factor 
Leaves-ground view factor 
F/^ w Leaves-wall view factor 

^ia_side«a Side of tree canopy above roof-sky view factor 
^ia_side«r Side of tree canopy above roof-roof view factor 


Ci Specific heat capacity of leaves J m~ 2 K 1 

Hi Sensible heat flux of leaves W m~ 2 

E/ Moisture flux on leaves kg m~ 2 s -1 

E root Transpiration from root zone kg m~ 2 s -1 

L„ Latent heat of vaporization J kg -1 

C w imm Specific heat capacity of 1-mm water depth J m~ 2 K 1 


then h c is set to fi[ c . The direct solar radiation absorbed 
by the leaves is 


5 


D l _ 
/ ~ 



(A7) 


/ground are the normalized shadow length, normalized 
building height, and normalized street-canyon width, 
respectively, as defined by Kusaka et al. (2001). The 
solar radiation absorbed by the leaves is a function of the 
leaf aspect ratio which is defined as 


and the direct solar radiation reaching the surface of the 
wall, ground, and grass/soil is 

S° l =(5f -5f*)kadow j (A8) 

■"height 

S% =/ rd (5f 1 - sf *) /ground ~ /shadow , and (A9) 

Aground 

Sf =f s ( S ^ ~ ) V°und~ Shadow ^ (A10) 

ground 

where / rd and f s define the fraction of the ground that is 
road and grass/soil, respectively, and /shadow, /height, and 


(T i = LAI*/), (All) 

where the effective leaf area index is 

LAP = 2.5 [1 - exp(— 0.4LAI)] , (A12) 

as described by Lee and Park (2008). 

The indirect solar radiation reaching the surface of the 
leaves if buildings cause no shade is defined by 

Sf* =S I c l f l [l-T(0,h / )], (A13) 

where S' is the downward indirect solar radiation at the 
top of the street canyon (i.e., roof level) (Lee and Park 
2008). 
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The indirect solar radiation flux absorbed by the 
leaves is (Lee and Park 2008) 


i-“/) 


S = 


T wg =1 -fl 


i-r( o,- 4 h b 


and 


(A14) 


= 1 —//[l — T(0,h {c )\, 


(A20) 

(A21) 


and the indirect solar radiation incident on the surfaces 
of the wall, road, and soil is (Lee and Park 2008) 


^=^ fl (S' l -sf), 

(A15) 


(A16) 

S's =f s F g „M l -S?), 

(A17) 


where F w „ a and F g ^ a are the sky view factors at the 
center of the wall and ground and /•)„„ is the sky view 
factor of the leaves at one-half of the height of the tree. 
The sky view factors are computed with the same algo- 
rithm as shown in Kusaka et al. (2001). 

The mean radiative transmissivities due to the tree 
canopy within the street canyon are calculated by (Lee 
and Park 2008) 


where T wa , r ww , T wg , and T ga represent the respective 
transmissivities between the wall and atmosphere 
above the canyon, between two building walls, be- 
tween wall and ground, and between ground and at- 
mosphere above the canyon. The UCM allows for 
solar radiation to be reflected twice within the urban 
canopy. The solar radiation absorbed by the building 
wall, road, soil/grass, and tree leaves is defined, re- 
spectively, as 

= (1 - <*„)[$ + -C + T wg « rd (S^ + Sfd )Vw 
+ T Wg“X +S f i )VH' 

+ J. (A22) 


, = 1 -// 

■ = Wz 


^~T[- A h b ,h ic 


i-n\h b ,\h b 



s r y 

(A18) 


(A19) 

and 


cl^ I nD^ 


C / 1 I c 


(A23) 


1 Wg-' .S' W W VP ' g*-*w J’ 

(A24) 


S] l =S ! ;+ s? 1 (1 - «,)[(! - r w > rd (^ + )F w „ g + (1 - r)a s (Sl l + Sf l )F t 


wg/ 
r/ I cD > 


w^g 


+ (1 - + (! - T wg )fr^J S C + S° l )F g „ w + (1 - T wg )f s a w (SC + S° l )F g „ w 

+ (1 - + S D J )F a „ w + (1 - r ga )« rd (4 + )F + (1 - r W(Sf + 1 )F a „ 1 , 


(A25) 


where subscripts w, g, and a represent the wall, ground, 
and atmosphere above the canopy, respectively, and F is 
the view factor, which is calculated in the UCM as given 
in Kusaka et al. (2001). 

Longwave radiation is emitted isotropically from the 
tree canopy. Some of the radiation emitted from the tree 
leaves reaches the surface of other leaves within the tree 
canopy while the remainder reaches the ground, build- 
ing walls, or atmosphere above the street canyon. The 
longwave radiation emitted from tree leaves that reaches 
the surfaces of other leaves, the atmosphere above the 


canyon, building walls, and the ground is defined as in 
Lee and Park (2008): 

L{ T =L ZT / z [l-r(0,A fc )], (A26) 

Lf t =0.5cr z (L zt -L{ T )F z „ a , (A27) 

W i 

L n =w°- 5a ^ L n - L n )( 2 “ F i»a ~ F i «*)• and 

An b 

(A28) 
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Lf T =0.5 CT/ (L /T -L' T ) J F ; „ g , (A29) 

where w is the width of the ground and L/j is defined as 

L, t =s l a l Tf, (A30) 

with e/ and F/ being the emissivity and temperature of 
the leaves, respectively. The UCM allows for longwave 
radiation to be reflected twice within the urban canopy. 


The net longwave radiation at the surfaces of the 
building walls, roads, grass/soil, and trees is obtained by 
incorporating the downwelling atmospheric longwave 
radiation; longwave radiation emitted by building walls 
and roofs, roads, grass/soil, and trees; and multiple re- 
flections in the urban street canyon. The net longwave 
radiation at the surface of the wall is 


L l l = LClF w~a + T wg(/id e rd F M.«* ff7 ri + f s s A 


^ ) + T. 


o f ryT 4 

ww w w^w w 


T w — 

*-l t al K 


+ ^wg[/r d(i - e ni )F w ~An +/,( 1 - r 1 + T ww(i “ e ^ F - 


W + AvgAvat/rdd - 


W' w vy 

+ /,(! “ es)F w „ g Fg~« LCl] + Wwgt/idC 1 “ e rd )F w „ g F w „ g e w crT* +/ s ( 1 - e s )F w „ g F w „ g e w vT w 


4 1 


+ Vugl/riC 1 - s w )F w „ w F w „ gBg aT 4 g + / s ( 1 - e w )F w „ w F w „ g e s (rT 4 ] + T w r wa ( 1 - e w )F, 


vy<->vy vy •<->« 


+ T T (1 — £ )i 7 
ww ww v vy/ v 


vy/ w vy vy<->vy vy 




L c l 

(A31) 


Here, L c 4 is the net longwave radiation entering the street canyon from aloft. The net longwave radiation at the 
surface of the road is 


=/rd e rdKa LCi + AvgV^W^ + L * - + T wg (l - sJF^L^ + r wg r wa (l - sJF^F^L^ 


/rd rd 
+ T, 


wg vy vy <->g 


wgv 


t F F 

wg ww g<->vy vy <- 


+ AvgAvgt^ 1 - e W ) F g«W F w«g e rd°' r r 4 d + /,(! “ 


F. K<r /;']}. 


g<->vy vy <->g 5 


(A32) 


The net longwave radiation at the surface of the grass/soil is 
ZU=f e { T L ei F +r eF crFl + Lf. - crF? + r (1 - e )F LX. + r r ( 1 -e )F F L ci 

s /s 5 ^ ga g<-># wg vy vy <->g vy /| ^ wg' vy/ g<->vy *T wg wa ' vy- 7 g<->w vy<->« 


wg 

+ t t F F s a 

wg ww g<-> vy vy«->vy vy 


T l + ^g^wgt/rdC 1 “ s w ) F g» w F w ^ g e ri <rT 4 li +/ s (l - eJF 




(A33) 


The net longwave radiation on the surface of leaves is 

rtf — j l r l r Z 4_ j ‘ 
*-/ ~ *-cl ^ *-12 ^ *- 


rd 2 


^2 F 


w 2 


-iy- 


(A34) 


where Lj. 2 represents the downwelling longwave building walls, respectively, that is absorbed by the 

radiation from the atmosphere above the street tree canopy. Quantities L l c2 , L l l2 , L[ A2 , L 1 s2 , and 

canyon that is absorbed by the tree canopy and take into account multiple reflections from the 

L l a , L l ti2 , L l s2 , and L l w2 represent the longwave radia- building walls, road, and soil in the urban street can- 

tion emitted from the tree canopy, roads, soil, and yon and are calculated as 


F'cl = ^{/ rd ^ V«( 1 - V) + fs LCl F g ~A ~ Tga) + d “ T ga )[r ga / rd (l 


T T f (1 — e )F F L ci l + (f-r )t (1-e )F 

ga Js' s' g*~>a a<->g J ' wa/ wa ' w' v 

+ f — b )F F L c1 ]+t (1 — t )(1— e )F 

J s'- s' vy<->g g*->a J wa v ww/' w' v 

+ f(l-8 )F F 

's'- w' g<->vy w< 


F 

w' w*-+a a*->w 


L C Ar(l-r)[f,(l-s Td )F K 


L cl 


F 

w*->w vy •< 


*A cl ]}, 


L ci +T (1 — r )[ f ,(1 — e )F F L cJ - 

wa'' wg^Ly r dv vy/ g-«-> vy vy 

(A35) 
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L\l = - v)Kd(i - e rd )V« L n +# 1 - e sK~ a L n] + (! - T wa)(l - 


+ 0- - zJK-A +/,(! - e S )^~/TT] + 0- T ww)( 1 - 

+ (i - -w g )Kd(i - s w ) F g ^ w L n +/,( 1 - OV» l "t]}’ 




r H 




(A3 6) 


-rd2 =^f[/id / g«»v c id a,:7 rf( 1 “ r wg) + /rd F g-a £ rd°' r r 4 d( 1 ~ + (1 “ T wa) T wg( 1 “ e w) F w^a F g^w e rd CTT rd 


F T wg(^ T ww^rd^ S w^ F w<r+w F w<-+g S rd CrF td T wg(^ T wg)^rd-^rd^ S w) F w<->g F g*->w e rd°"^rd 

r 4 


+ T wg( 1 “ ’ , w g )// rd ( 1 - e w )^«H- / M-«g B i d 0 ' :r r d ]. 


(A37) 


L *2 = “ T wg) + 4 F *«a 8 *°' 7 ?( 1 “ T ga) + C 1 “ T wa) T wg( 1 “ 


+ T w g (! “ at 1 “ 8 iv ) F w ** w F w <^e B s CrF s + t we (1 - T w M(1 ~ eJF w 


w' w<-+w w*-+g S s Wg v 


Wg 'J S J S' 




+ T wg(! “ T wg )//,(! “ B J F g»w F w»g B s (T Ts l and 


(A38) 


L 'wl =^{f I d F g~w B w< rT l( 1 ~ r wg) + / I /7 g« w e w °' 7 2( 1 - T „ g ) + F w „ w e w <rl$( 1 - T ww) + - T wa) 


+ 0 - T ga)[ T wg/rd( 1 - s rd ) F g^ a F w ^g B w (lT i + ^wg/4 1 “ e s ) F g »a F W ~g B w‘ xT $>~\ 


+ n — t It n - £ YF F 

v wa 7 ww^ w 7 u><->a w 

p4i 


+ f (1 - 8 )F F 8 0-711 + T (1 - T )(1 - 8 )F 

J S' S' g<->w W^g W Wi WW' WW 7 ' W' l 


T w + T wg(! - T wg)[/rd( 1 “ E rd ) F g ~ w F w ~g B w ° T w 

4 


WW' 

r 4 


w w*^w F w*-^w B w t ' F w 
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If trees are taller than the adjacent buildings, then the 
heat and moisture fluxes are calculated for the portion of 
the tree canopy above the roof height and tree shading 
of the roof is considered. The direct solar radiation and 
indirect solar radiation that reaches the leaves above the 
roof height are 


absorbed by the tree canopy entering the top of the tree 
canopy are defined as 

Sg* = 5 ‘ a ^ ^ and (A42) 


and (A40) 

€ =S'aW-mh !a )], (A41) 

where SP * , 5?* , Sf* , S F , and /if a are the downward di- 
rect solar radiation that reaches the surfaces of the 
leaves, downward direct solar radiation from the atmo- 
sphere above, downward indirect solar radiation 
reaching the surfaces of the leaves, downward indirect 
solar radiation from the atmosphere, and height be- 
tween the treetop and building roof level, respectively. 
The amounts of direct and indirect solar radiation 


Sfa = >S|a ^ a ‘\ (A43) 

°7 

Solar radiation is also transmitted through 
the sides of the tree canopy above the roof. The 
respective direct solar radiation and indirect 
solar radiation absorbed through the side of the tree 
canopy are 


cD ^ 

^la_side 


^shadow roof S^ a (1 OL^) 


and 


(A44) 
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S /l = 

la _side 




(A45) 


where / shad0 w_roof and F laside „ a are the normalized 
shadow on the roof and the view factor between the side 
of the tree canopy and the atmosphere above, re- 
spectively. The amount of solar radiation reflected from 
the roof and absorbed by the tree canopy is 


where Sj 1 , a r , and F laside „ r are the solar radiation ab- 
sorbed by the roof, albedo of the roof, and view factor 
between the roof and the side of the tree canopy, re- 
spectively. So, the solar radiative flux on the tree canopy 
above the rooftop is 


r r/T I cF i rifi 

^la - ‘-’la ' 1 ~ ‘ha ^ ‘ha.side 


■ 

^la_side 


■ s ^ 

°la_r- 


(A47) 


s J- = 

^la-r 


S}l[l-T(0,h u )]f l a r Q.-a l )F l 


la_side<-»r 


(1 - a r )a l 


(A46) 


The solar radiative flux on the roof is 


SP = 


^roof ^shadow 


(1 -a r )S? + {l- 


roof 


^roof ^shadow 


roof 




)]} 


+ 5^(1 - a r )F r ^ a + S'*(l - «,.)(! - F_ fl ){ 1 -/ z [ 1 - 1(0,/.,)]}, 


(A48) 


where h a ^ rdkca is the distance of the path of the direct 
radiation through the tree canopy from the center of 
the shadow on the roof (/ s hadow ro „[/2). The first term in 
Eq. (A48) is the amount of direct solar radiation that 
reaches the roof with no obstructions and is absorbed by 
the roof, the second term represents the transmitted 
direct solar radiation through the tree absorbed by the 
roof, the third term is the indirect solar radiation that 
reaches the roof with no obstructions and is absorbed by 
the roof, and the fourth term is the indirect solar radi- 
ation transmitted through the tree canopy that is ab- 
sorbed by the roof. 

When the trees are taller than the building height, the 
amount of longwave radiation entering the street can- 
yon is 

L ct = L ai {1 _ /;[1 _ r(0,fi fa )]} + Lf at , (A49) 

where LjU is the amount of longwave radiation emitted 
from the tree canopy above the rooftop that enters the 
street canyon, defined by 

Z4 t =0.5cr ; (L laT -L£ T ), (A50) 

where L la j is emitted from the tree canopy above the 
roof defined by 

L la T = e / CT X ( A51 ) 

and /.[:], is the amount of longwave radiation that is 
emitted from the tree that reaches the surface of other 
tree leaves above the roof height: 


^!at= L laT //[l- 7 ’(0,/i fa )], (A52) 

where F la is the temperature of the leaves above the 
roof. 

The longwave radiation flux when the trees are taller 
than the adjacent buildings is defined as 

Lj 1 = * r (L^F r „ a + L a * (1 -F r „ a ){l-f,[l-T(0,h {a )]} 

-a7? + 0.25^(L -Li* t ) 

ia 

x (2 _ ^ir ~ -Fi r „ r )-F r „i r _ s i de )> (A53) 

where the first term is the amount of longwave radiation 
from the atmosphere that is not obstructed by the tree 
canopy and is absorbed by the roof, the second term is 
the amount of longwave radiation from aloft that is 
transmitted through the tree canopy and absorbed by 
the roof, the third term is the amount of longwave ra- 
diation emitted by the roof, and the fourth term is the 
amount of radiation that is emitted by the tree canopy 
and absorbed by the roof. 

The longwave radiation flux in the tree canopy 
above the rooftop is the sum of the amount of radia- 
tion that is emitted from the road, grass/soil, building 
walls, building roof, tree canopy in and above the 
canyon, and atmosphere above the trees that is cap- 
tured by the tree leaves minus the amount of radiation 
emitted from the tree canopy above rooftops, which is 
defined as 
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